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We present a hidden abelian extension of the Standard Model including a complex scalar as a dark 
matter candidate and a light scalar acting as a long range force carrier between dark matter particles. 
The Sommerfeld enhanced annihilation cross-section of the dark matter explains the observed cosmic 
ray excesses. The light scalar field also gives rise to potentially large cross-sections of dark matter 
on the nucleon, therefore providing an interesting way to probe this model simultaneously at direct 
and indirect dark matter search experiments. We constrain the parameter space of the model by 
taking into account CDMS-II exclusion limit as well as PAMELA and FermiLAT data. 



I. INTRODUCTION 

The existence of dark matter (DM) in the present Universe has been firmly supported by a range of evidences [Ij . The 
prime among them are galaxy rotation curves, large scale structure, cosmic microwave background and gravitational 
lensing. However, the identity of DM within the Standard Model (SM) of particle physics is still missing. Therefore, 
its experimental verification is expected to be a new discovery and a strong indication for physics beyond the SM. 

Recently the Cryogenic DM Search (CDMS) Collaboration in the Soudan mine reported the observation of two 
events compatible with a positive dark matter detection at 1.64 a confidence level (C.L.). The DAMA/LIBRA [3] 
experiment claims an evidence of DM in its modulated signal at 8.8 a C.L. Several other direct detection experiments 
are running and setting upper bounds, including EDELWEISS-II g], ZEPLIN-II (5|, XENON-10 [6] and Cresst-III [7]. 
The Xenon collaboration will soon release the data from the first run of the XenonlOO experiment jS] and superCDMS 
is planned. 

A significant amount of effort has also been devoted to detection of DM through indirect searches. For example, 
the PAMELA collaboration [9_ reported an unexpected rise of the positron fraction compared to that of the galactic 
background at energies above 10 GeV, while confirming the earlier results of AMS [TU] and HEAT [TT|. Similarly the 
HESS [T2] and Fermi Large Area Telescope (FermiLAT) [13J collaborations also reported an excess of electron plus 
positron flux with respect to the galactic background at energies above 100 GeV, but without confirming the spectral 
features observed by the balloon-based experiments ATIC [14] and PPB-BETS [15] . It has been widely interpreted 
that DM could be a viable candidate for the observed cosmic ray anomalies, although they could be explained by 
astrophysical sources |16| . 

In light of the above experimental results, several models have been considered in the literature, which leave 
signatures at either direct[T7] or indirect |18l ITU] DM searches. Typically, for a given model, the predictions for 
direct and indirect signatures of DM depend on different parts of the parameter-space, and the derived constraints 
thus do not overlap. However, in some models the same couplings are responsible for both the scattering of DM on 
the nucleon and large annihilation cross sections, in which case an interesting complementarity between direct and 
indirect searches exists |20l |2"T] . 

A simple possibility is to consider singlet extensions of the SM, in which the DM is a singlet scalar [22 , coupled 
to the SM Higgs particle through the so-called Higgs portal [23 , namely the Higgs to DM quartic coupling (see also 
doublet extensions of the SM [24]) ■ Depending on the strength of this portal, the singlet can account for the observed 
relic density, Vl^uh 2 ~ 0.1 [25J. Furthermore, through the Higgs portal coupling, DM scatters with the nucleon and 
is thereby constrained by direct searches, as well as annihilates into SM fermions which can be observed at indirect 
detection experiments [201 12T] . 

The anomalous positron and electron fluxes observed by PAMELA and FermiLAT require a large enhancement 
of current DM annihilations. In their minimal versions [25], the singlet extensions cannot reproduce such features. 
Furthermore, in these models, the stability of DM is supported by an ad hoc discrete symmetry. In this paper, we 
study a model which naturally solves these issues. We introduce a hidden sector gauged under an Abelian U{1)h [2"B"] - 
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|27| . containing two complex scalars $ and S. While all SM fields are hidden sector singlet, the extra scalars are 
singlet under the SM but charged under U(1)h- This model provides all the ingredients for a viable DM model with 
potentially large direct and indirect detection signals. 

The paper is organized as follows. In Sec. [IT] we present the model, with particular attention to the mass spectrum 
and mixings in the scalar sector. The model parameter space is then examined in Se c.|III| and constrained by requiring 
the relic density of the DM candidate to be in the WMAP7 range. Sections |IV| and | V| describe the phenomenology of 
the model in light of the present DM searches: CDMS-II, XENON10, PAMELA and FemiLAT. First, in Sec. [TV} we 
discuss the direct detection bounds constraining the hidden sector parameter space and the interplay with the indirect 
detection. Then, in Sec.|Vj we investigate the indirect detection bounds, and give the results in terms of positron and 
electron excesses. The conclusions are presented in Sec. |VI| 

II. THE MODEL 

Assuming a non-trivial charge assignment under the hidden Abelian U(1)h for the extra scalars S and $, qn(S) = 3 
and <7h ($) = 2, the most general scalar potential is given by: 

+f HS H^HS^S + fB+rfH&Q + fsi&S&Q . (1) 
The hidden sector couples to the SM via Higgs portals, as schematically depicted in Fig.Q. We assume that all 
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FIG. 1: Schematic representation of Higgs portal couplings to a hidden sector. 

parameters are real and positive. The field <i> acquires a nonzero vacuum expectation value (vev), which triggers the 
breaking of U(\)h to a remnant Z 2 symmetry under which all fields are even but S. The latter does not develop any 
nonzero vev and hence can be a dark matter candidate, being stabilized by the Z 2 symmetry. 

It is remarkable that almost all the parameters in this model are well constrained by both direct and indirect 
searches, as we will see in great details below. The relic density of S is obtained through its annihilations into H 
and $ via fns an d fs</>> respectively. These two couplings also enter in the spin-independent cross section of DM on 
nucleon and hence are strongly constrained by direct DM searches. The field <i> gives rise to a long range attractive 
force between S particles, thus boosting the current S annihilations, while keeping the relic abundance unchanged. 
As a result fs$ gets strongly constrained from indirect DM searches. After $ acquires a vev, it mixes with the SM 
Higgs through fu<t>- 3> is destabilized and consequently decays into SM fermions through the same coupling. While 
fjj<t> should not be too small for $ to decay before the onset of big bang nucleosynthesis (BBN), a too large H-<& 
mixing is excluded by constraints coming from LEP-II on the ratio of the invisible-to- visible Higgs decay |28| . 

A. Masses and Mixings of Hidden Sector Fields 

From the scalar potential in Eq. (JTJ , let us derive the quantities relevant for this study. First of all, vacuum stability 
requires, besides positive quartic couplings, that A\}j\g\^-\- fHsfH4>fs<p > + ^(pfns + ^sfH<t>- The electroweak 

symmetry breaking occurs when the SM Higgs acquires a vev (H) — v/y/2, while U(1)h is broken to a surviving Z 2 
symmetry when $ acquires a vev ($) = u/y/2. In the unitary gauge, the quantum fluctuations around the minimum 
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are parametrized as 

Hss ^W +h ) and $ -7i (u + 0) ' (2) 

where h and <p are physical real scalars, the unphysical degrees of freedom being eaten by the longitudinal component 
of the SM gauge bosons and of the Z' associated with U(1)h- Minimization of the scalar potential in Eq. ([I]) enforces 



v _ f 2 Mh A</> ~ m| fH4> \ u _ f 2^4 Ah - n 2 H fu4> \ 



(3) 



This minimum is the global one if fjj, < 4AhA<£. The two real scalars /i and <j> mix with each other and the mass 
matrix in the basis spanned by (h, <fi) is given by: 

M 2 =(f Hv2 f J!+™) . (4) 

Assuming /h^ <C 1, the mixing angle between h and is suppressed: 

e "~- 2<X f yi ( 5 ) 
2 (A_hV - A^u ) 

In terms of (h, 0), the mass eigenstates h\ and /12 read: 

h\ ~ h + 9 rn <fi and fi2 ~ — # m ft • (6) 

Consequently /ii is mostly the SM Higgs field, while /12 is the light scalar. Their respective masses are: 

M 2 hi ~ 2A H « 2 and Af 2 2 ~ 2A 0U 2 ^1 - . (7) 

Because of the small mixing, the current experimental bounds [281 129] on the SM Higgs mass apply on Mft x . Hereafter, 
in all numerical evaluations we take Mh 1 = 120 GeV. However, we note that our results are quite insensitive to the 
Higgs Mass. The light scalar field h,2 is the main product of S annihilations at the present epoch. In order to avoid 
overproduction of high-energy gamma rays from the decay of h 2 , M^ 2 < 2m^ is required: we take M^ 2 = 240 MeV 
as numerical reference value. We then have A#u 2 3> A^u 2 . 
The DM mass is given by: 

Ml =4+ ^ + >|V , (8) 

which is varied from 10 GeV to 1 TeV in the following. 

The hidden sector is gauged under the Abelian U(1)h and contains an extra gauge boson Z' . Although SM particles 
are singlet under U(1)h, through the kinetic mixing of Z' with the hypercharge gauge boson, they can couple to Z' . 
The relevant part of the Lagrangian then reads: 

Lz< = -\F«F H ^ - \F* V F H ^ + - 2ig H Z'^\ 2 + - UgnZ'^ , (9) 

where gn is the hidden sector coupling constant, FX, and F^ v are the field strength tensors associated with U(l)y 
and U(1)h respectively and x parameterizes the kinetic mixing between the U(l) symmetries. 

For an invisible Z' , the kinetic mixing is expected to be x ~ 10~ 4 -10~ 2 at the electroweak scale. The mass of Z' is 
given by Mz> = 2gnu and we take 600 GeV as a reference numerical value. The parameter u is then lower bounded 
by the requirement of perturbative U(1)h gauge coupling: u > 85 GeV. For the chosen value of the Z' mass, the 
electroweak precision measurement constraints the kinetic mixing to be x < 0.036 [50]. Through Z-Z' mixing, the 
DM can scatter on nuclei. The corresponding cross section may exceed the current experimental bounds for x ^ 0.01. 
When Mz> < Mg the dark gauge boson is produced in S annihilations. For large x> this would lead to significant 
antiproton fluxes in the cosmic ray which is in contradiction with PAMELA data |31| . In the following we assume 
X ~ 10~ 4 and so Z' effects are negligible. 
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Assuming Mz> = 600 GcV with gn ~ 0{1) while Mf l2 = 240 MeV implies a strong tuning of the model parameters, 
with /x^ = 120 MeV <C « ~ 300 GeV and ~ 10 -7 . Reducing this tuning can be achieved by either increasing 
the hidden scalar boson mass, or by lowering Mz' as well as Ms- However, in those cases, predictions of the model 
turn out to be very different from those we discuss hereafter. By increasing Mh 2 we possibly face an antiproton 
overproduction in DM current annihilations, in which case the corresponding cross section should be suppressed 
compared to PAMELA requirements. In this case, still, the model can afford a viable DM candidate but lacks 
predictivity. Another way to reduce the tuning is to lower Mz> or Mg. With a lighter DM, the PAMELA cosmic ray 
spectrum cannot be accounted for, while with a lighter Z' the model is in conflict with direct detection experiments. 
Therefore, in the following, we assume a certain amount of tuning, allowing the model to be predictive and probed 
by both direct and indirect detection experiments. 



B. Astrophysical and collider constraints 



We now discuss some relevant constraints on the hidden sector coming from astrophysics and the electroweak pre- 
cision measurements. The scalar field hi decays to SM fermions through <$>-H mixing; its decay rate is approximately 
given by: 

r*,* £ (10) 

2m f <M h2 

The lifetime of hi is then estimated as 

Tu. ~ 0.1 s x 

M h2 \10 

Thus, demanding that h% decays before the onset of BBN 77^ ^ tbbn ~ 0.1s [32] , 9 m should be bigger than 10~ 7 
Equivalently a lower bound on the Higgs portal coupling fjj^ is inferred: 



T,„~0.1»X^X(-^1 . (II 



Assuming that hi dominantly decays to fi + fi~ , a strong constraint on 6 m can be obtained from the £>-meson decay. 
From the branching ratio BR(B — > /i + fi~ X) < 3.2 x 10~ 4 [27 ,[33J, one infers an upper bound on 8 m to be \6 m \ < 
10~ 2 P5- This in turn gives 

fe<10-xf^V^V. (13) 



V 120 GeV 

Since /12 is very light, the mixing angle never saturates the above bound. For the same reason, the electroweak 
precision measurements on the S, T and U parameters do not receive any significant contributions from hi |35| . 



III. ANNIHILATION CROSS-SECTION OF THE 5" DM AT FREEZE-OUT AND AT PRESENT EPOCH 



A. Relic density 



In the early universe the S particles maintain their thermal equilibrium through the scattering processes: S'S — > 
hihi and S^S — > SM particles. In Fig. we display ratios of the dominant annihilation channels with respect 
to the total annihilation cross-section <7tot> as a function of the DM mass. In case of fus >> fs<p (left panel) the 
dominant channel is S — > // (blue dashed curve) until the W + W~ (black dot-long dashed solid) and ZZ (black 
dot-dashed) thresholds are reached. These two processes dominate the behavior of S annihilations up to 120 GeV, 
when S — >• h\h\ (green dashed curves) takes the upper hand. The situation is reversed as soon as fg^ > fjjs 
(central and right panels). The annihilation into h\h\ is competitive only for equal Higgs portal couplings, while for 
fscj> — 10 Jhs, S^S — >• hihi (red solid curves) dominates in all the DM mass range, apart from the Higgs pole. 

As the universe expands, the temperature of the thermal bath gradually falls; at Tjj ~ Ms/25, S gets decoupled 
and starts redshifting. The relic abundance of DM can be evaluated by solving the Boltzmann equation for the S 
number density: 

^ + 3n s H(t) = -(a\v lel \) (n| - n c f) , (14) 



5 




fHS = 0-1 fs« 







-- 












% : 



M s (GeV) 



M s (GeV) 



FIG. 2: Ratios of the leading annihilation channels of S on the total annihilation cross-section are shown as a function of Ms: 
in red solid (S*S -> h 2 h 2 ), in green dashed (S*S -> hihi), in blue dotted (S f S -> //), in black dot-dashed (S^S -» W+W~) 
and in black dot-long dashed (S'S — > ZZ). 



where H(t) is the Hubble expansion parameter and 
section. 



'<7|i> re i|) is the thermal average of the total annihilation cross- 



The model has been 
The parameter space is chosen as follows. The physical Higgs 



To accurately evaluate the DM relic density, we solve Eq.(|14| by using micrOMEGAs 
implemented in micrOMEGAs through FeynRulcs [37_. 
masses are fixed at = 120 GeV and M^ 2 — 240 MeV, while the other relevant parameters are varying randomly 
within their allowed ranges: u from 10 GeV to 1 TeV and fi s from 1 GeV to 500 GeV. The portal couplings fs$ and 
fns vary from 10~ 3 to y/4n, the perturbative upper bound, while fg^ varies from 10~ 7 to 10~ 4 . We also constrain 
the lifetime of h 2 to be less than 0.1s. 

Demanding that the relic abundance of DM should satisfy the WMAP7 constraints [25] at 3a C.L., given by 



n 



DM' 



0.0941 -0.1277, 



(15) 



we show in Fig. ^ the resulting scatter plots in the plane of fs<j> and fjjs against Ms- 

The requirement of having the correct relic density fixes the balance between the two Higgs portal couplings fs<f, 
and fas and hence between the different annihilation channels. In the left panel of Fig. as long as M§ < M^/2, 
we clearly see that the S — > h 2 h 2 channel dominates. Near the Higgs resonance, Ms ~ M^/2, lower values of the 
portal couplings are allowed, but at the pole, annihilations are too efficient and the relic density gets suppressed. For 
< Ms < Mw, the DM dominantly annihilates into light SM fermions, mostly bb pairs. For Mw < Ms < Mh., 
the S^S — > W + W~ and S —> ZZ channels, mediated by hi, strongly constrain fns, as can be seen in Fig. pj) 
(right panel). For larger DM masses, the annihilation channels S^S — > h\h\ and Sr S — » it are also allowed. For such 
large DM masses, fns is more constrained than fg^. The latter can indeed take values up to the perturbative bound, 
while the former is upper bounded at 1. 

Resuming, on the whole range of DM mass, /s^ typically takes values of the same order or larger than fns- This 
corresponds to middle and right plots of Fig.Q, and implies that the DM is preferably annihilating into h 2 , as can 
be inferred from the approximate cross-sections: 



(a\v lel \)(&S^h 2 h 2 ) 
(a\v lel \)(&S^hihi) 



f 2 



64ttMJ 



f 2 



HS 



64tt M| 



Ml 
Ml 



1/2 



(16) 



Thus we see that both portal couplings fns and fs<p are constrained from the requirement of obtaining the right 
DM relic density. Interestingly, these couplings also enter in the t-channel Higgs-mediated scattering Sf — ► Sf, which 
is relevant for the direct detection of DM, as we discuss in the next section. 



B. Sommerfeld enhancement 



The observed relic density of DM can be obtained for an annihilation cross-section (cr|?; rc i|) D ~ 3 x 10 _26 cm 3 s _1 
at the decoupling epoch, To ~ Ms/25. However, at the present epoch, such a cross-section is too small to explain 
the anomalous cosmic ray fluxes observed by PAMELA and FermiLAT. Several mechanisms have been invoked in 
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FIG. 3: Scatter plot for the relic abundance of S. In the left (right) panel the solutions that satisfy Odm/i 2 = 0.0941 — 0.1277 are 
shown as black points in the plane Ms versus fs4> (fas)- The green (gray) region denotes points with boost factors compatible 
with indirect DM searches. 



the literature to boost it. A natural enhancement may occur due to the local overdensities of DM in clumps |38| . 
but estimates show that the resulting boost of the annihilation cross-section is too small to account for the observed 
cosmic ray fluxes. Two different mechanisms, arising from particle physics, also exist: either a Breit-Wigner resonance 
of the annihilation cross-section [27 ,[39 , or the Sommerfeld effect [3D]. I n this model, the latter naturally occurs 
because of the presence of h^- The light scalar h% acts as a long range attractive force carrier between the DM 
particles. If S kinetic energy is small enough, the attractive interaction becomes relevant and induces an enhancement 
of the annihilation cross section. Denning the reduced DM bound-state wave-function as tp, the corresponding boost 
is computed by solving the radial Schrodinger equation: 

ij}(r)" - M s V(r)ip(r) + Mj f3 2 t/j(r) = 0, (17) 

where f3 = v re i/c is the DM relative velocity, and V(r) is the attractive Yukawa potential: 

V(r) = --e- M ^ r , with a=(% U ) 2 ~ . (18) 



r 



V M s J 8tt 



We solve Eq.(17l using the boundary condition ij/(oo) — \Msf3 tp(oo). The Sommerfeld boost is then given by 

■0(oo) 



S< 



(19) 



To consider the enhancement at present time, one should integrate over the velocity distribution of DM in the Earth's 
neighborhood: 

(S e )= AUm / dv Tcl -f e-^/*° S e (tw ,a,M h2 /M s ) , (20) 
Jo v 

with mean velocity vq — 220 km/s, escape velocity v csc — 650 km/s and a normalization factor A/" n0 rm for a smooth 
maxwellian halo. The boost factor is then only a function of u x fs$, Ms and M^ 2 . Therefore, for a given DM 
mass, (S e ) is degenerate with respect to fss and u, the hidden sector breaking scale. The effective DM annihilation 
cross-section then reads 

(a\v rcl \) = (S e ) x (a\v Iel \) D . (21) 

In Fig. (|) the region of the parameters giving rise to boost factors from 20 up to ^2000 is displayed by a green 
(gray) region. This is the range of enhancement required to explain the observed cosmic ray anomalies for a DM of 
mass (100 - 1000) GeV, compatible with current cosmological and astrophysical constraints, as discussed in Sec. [V 
On the left panel we see that only large values for fscj, are allowed, ranging from 10 _1 up to This is expectec 

by the behavior of the Sommerfeld enhancement oc a//3. As mentioned before, the requirement of having a good relic 
density fixes the balance between the two portal couplings. In the right panel of Fig. the values of fns in the 
green (gray) region correspond to the values of fs<f, giving rise to a correct boost factor, although the Sommerfeld 
effect does not depend on fas- For 100 GeV < M5 < 200 GeV, the largest values of Jhs are n °t allowed because 
they dominate S annihilations, which in turn corresponds to lower values of fstj>- 
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FIG. 4: Scatter plot for the direct detection of S particles in the plane Ms versus asi- Solutions with a relic density in the 
WMAP7 range are depicted by black points; the red region (dark grey) is excluded at 90% C.L. by CDMS-II, while the blue 
region (light gray) is excluded at 90% C.L. by XENON10. The solid blue curve is the forecast upper bound for the XenonlOO 
experiment. The green region (gray) denotes points that in addition provide a correct boost factor. 



IV. CDMS-II EVENTS AND S DARK MATTER 



A. Experimental upper bounds 

The CDMS collaboration has recently published the analysis of the final run of the CDMS-II experiment After 
the background subtraction and cuts, two events survive, respectively at 12.3 keV and 15.5 keV recoil energies. The 
significance of the two events being a DM signal is at 1.64cr C.L., namely there is 23% probability that these are 
of more common origin, such as cosmogenic or neutron background. In Ref. |41| . it has been shown that if the two 
events are taken to be DM signal, the 1.64ct region will prefer light dark matter candidates, with an upper bound 
on the mass around 60-80 GeV. If we consider the 90% C.L. region to constraint the DM mass, then only an upper 
bound can be set. We approach the analysis of the two events of the CDMS-II in a conservative way, finding a 90% 
C.L. exclusion limit with the maximum gap method |42j . 

In addition the upper bound from the XENON10 experiment |6J is taken into account, which is the most constraining 
one together with CDMS-II, in the case of spin-independent elastic interaction on nucleon. Using the data of the 2007 
run, reanalyzed as in Ref. |43| , and the maximum gap method, we infer the 90% C.L. upper bound. We also consider 
the prediction for the first run of the XENON100 experiment, following Ref. |S]. 



B. The elastic spin-independent S cross-section 

The interaction of S on the nucleon gives rise to a coherent spin-independent elastic scattering, mediated at tree 
level by the scalars of the model, the SM Higgs particle and hi. This cross section reads: 



Mn _2 f 2 1 fHS 1 fs<pU( 

2 Ml 2 Ml 2 . 



(22) 



where m n is the nucleon mass and fi n is the nucleon-DM reduced mass. The parameter /„ indicates the effective Higgs 
nucleon interaction, f n m n = (n\ Y] m q qq\n), where the sum runs over all the quark flavor. The /„ factor introduces 
hadronic uncertainties in the elastic cross section: its value vary within a wide range 0.14 < /„ < 0.66, as quoted in 
Refs. [33]. Hereafter we take f n to be 1/3, the central value. 

From the Lagrangian in Eq. ^ , the S particle couples directly to the hidden gauge boson and through the kinetic 
mixing it communicates to the SM fermions. This results in an additional DM-nucleon scattering. Neglecting terms 
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FIG. 5: Left panel: maximal ctsi value in the plane fus and u x f$cj> for Ms = 400 GeV, 6 m = 1-9 x 10 6 (blue dashed curve), 
9m = 5 x 10 -6 (dotted orange line) and 9 m = 10" 5 ( green long dashed line). The region plane on the right-hand side of each 
curve is excluded by CDMS-II at 90% C.L.. Right panel: scatter plot in the plane of u x fs<p vs Ms- Black points give the 
required relic density, while the points in the green (gray) region give rise to a boost factor compatible with PAMELA anomaly. 



of order % 3 , the corresponding spin-independent cross section is approximately: 

af{Sn ^Z'^ Sn) * % ^ ' X 2 sin 2 9 W (A + 2Zf , (23) 



where A and Z are the mass and the atomic number of the nucleus and g is the SU(2)^ coupling. The cross section is 
slightly dependent on the DM mass, through the reduced mass /i„. For Mz' = 600 GeV, the values of u we scanned 
over and for a maximal kinetic mixing \ — 0.036, the cross-section value ranges from 10 _45 cm 2 up to 10 _43 cm 2 . As 
stated above, we fixed \ = 10~ 4 , and consequently the cross section varies from 10 _50 cm 2 up to 10 _48 cm 2 , and is 
therefore a negligible correction to hi and hi contributions. It is then possible to constrain Higgs portal couplings 
thanks to direct detection searches. 

The behavior of ctsi as a function of the Higgs portal couplings is rather involved. Owing to the large mass differences 
between the two scalars, c/. Eq. ([5]), the mixing angle between hi and h\ is small, of the order of 9 m ~ 0(1O~ 6 -1O~ 5 ). 



Even though 9 m <C 1, the second term in Eq. (22 1 is not negligible in the whole S mass range. With respect to the 
standard Higgs exchange, the hi contribution is enhanced due to its small mass, which compensates the smallness of 
its coupling to the nucleon. Note that this result is valid in general in models where a scalar with mass lighter than 
1 GeV, typically the light force carrier of the Sommerfeld enhancement mechanism, mixes with the SM Higgs |20|[2T| . 

The predictions for ergi are shown in Fig. Q as a function of Ms, with all the points having a relic abundance in 
the WMAP7 range. The cross-section is enhanced respect to the standard Higgs exchange: a large region (red region) 
of the parameter space is not allowed by CDMS-II experiment. In the low mass range a portion of the parameter 
space is incompatible with the XENON10 upper bound, denoted by the blue (light gray) band. The green (gray) 
region describes the parameter space that leads to a large sommerfeld effect and boosts (cr|?; rc i|)£>. As described 
in the previous section, a large Sommerfeld enhancement calls for large fscj, coupling. We can therefore constraint 
the parameter space yielding such a large boost factor with the direct detection bounds. As shown in Fig. Q, a 
large portion of the green (gray) region is excluded, but nonetheless a large portion is found compatible with direct 
detection constraints. We also show the prediction for the upcoming XENON100 first run (blue line): while it can 
probe a bigger portion of the hidden sector parameter space compatible with indirect detection, a large part can give 
DM to nucleon cross-section below the expected sensitivity. 

On the left of Fig. ((5}, we show an illustrative example of the balance between the two contributions in Eq. (1221) , 
for a fixed DM mass Ms = 400 GeV. The maximum value of the elastic cross section compatible with CDMS-II 
is plotted, as a function of fns an d u x fs^. In this plot, three different values of the mixing angle are depicted, 
9 m = 1.9 x 10~ 6 (blue dashed line), 5 x 10~ 6 (orange dotted line), and 10 -5 (green long dashed line). First of all, for 
Jhs ~ 1) the SM Higgs exchange dominates and ogi is rather insensitive to 9 m . As fjjs decreases towards smaller 



values, the hi contribution becomes the leading one. For fus *C 1, from Eq. (22 1 we see that <tsi omv depends on 
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the product 9 m x u x fg^. For fns ~ 10~ 2 , we then have that different values of the mixing angle imply different 
maximum values of u x fs<f,. Notice that, from Eqs. (fl8| and (20 1, the Sommerfeld enhancement is sensitive to the 
same combination u x fg$. The right panel of Fig. T5J) shows - the points giving the right relic abundance in the 
plane u x fg^ vs Mg, where again a green (gray) region highlights the values of interest for the indirect detection. 
For fjjg = 10~ 2 , a small mixing angle 8 m = 1.9 x 10~ 6 allows for u x fg$ ~ 410 GeV, which is compatible with 
indirect detection constraints. Indeed, we see that for Mg = 400 GeV, u x fg^ should be in the 300 — 500 GeV range. 
Increasing 9 m calls for smaller values of u x fg^, and consequently lower (S e ) can be obtained. For 6 m = 5 x 10 -6 , 
we see that u x fg$ < 160 GeV, a value for which the boost factor is too small to account for the whole cosmic ray 
excesses. The situation is even worst for 6 m = 10~ 5 for which u x fg^ < 80 GeV. 



V. ELECTRON, POSITRON AND ANTIPROTON FLUXES FROM S f S ANNIHILATION 



As described in Sec. Ill the annihilation of S — > h 2 h 2 and S 1 * S — >■ h-Ji^ will generate both positrons and 



antiprotons in the present Universe and could have been detected by various experiments such as PAMELA, HESS 
and FermiLAT. In the case S annihilates into h 2 , the final products are dominantly muons and antimuons resulting 
from the decay of h 2 , while in the latter case the final products are mostly hadrons. From Eq.(16), the ratio of S 
annihilation cross sections into hi to h 2 is 



(a\v I e l \){&S-+h 1 h 1 ) ^ ffHS 

(a\v lel \)(StS^h 2 h 2 ) ~ [fa 



(24) 



As we see from Fig.(|3j), we typically have fg^ >, 10 fjjs when we require that the annihilation cross-section is enhanced 
by a non-negligible boost factor. This is sufficient to suppress the antiproton flux over the positron one. Therefore, 
in what follows we will focus on the production and propagation of positrons in the Galactic medium. 



A. Production and propagation of positrons 

From S annihilations, h 2 is produced which then decays to muon and antimuon. They ultimately decay to electrons, 
positrons and neutrinos. As a result, equal numbers of electrons and positrons are produced from the annihilation of 
S particles. However, the background flux of electrons in the Galactic medium is significantly larger than the positron 
one. Therefore, it is easier to find signature of DM, if any, in the Galactic positron flux. 

Once the positrons are produced in the Galactic halo where the DM concentration is large, they travel under the 
influence of the Galactic magnetic field which is assumed to be of the order of a few microgauss. The motion of 
positrons can then be thought of as a random walk. In the vicinity of the Solar System, the positron flux can be 
obtained by solving the diffusion equation [45 

-f e+ (E,r,t) = Q(E,r) + K e+ (E)V 2 f e+ (E,r,t) + -[b(E)f e+ (E,r,t)} 1 (25) 

where f e + (E, r) is the number density of positrons per unit energy, E is the energy of positron, K e + (E) is the diffusion 
constant, b(E) is the energy-loss rate and Q(E,r) is the positron source term. The latter, due to S annihilations, is 
given by: 

Q(E,f) =n%[?){a\v iel \) dI ^. (26) 

In the above equation the fragmentation function dN e + jdE represents the number of positrons with energy E which 
are produced from the annihilation of S particles. We assume that the positrons are in steady state, i.e. df e + jdt = 0. 



Then from Eq. ( 25 1 , the positron flux in the vicinity of the solar system can be obtained in a semi-analytical 
form @5M] 



Ms dN f 



$e+(S ' f0) = ^bW) {ns)2&{aM) J E dE'^I(X D (E,E')), (27) 



where \o{E,E') is the diffusion length from energy E' to energy E and I(Xd{E, E')) is the halo function which is 
independent of particle physics. An analogous solution for the electron flux can also be obtained. 
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FIG. 6: Positron fraction from S 
point 1 (2) quoted in Table [I] 



2/x + /j, at M$ = 400 GeV. The red solid line (black dashed line) stands for the illustrative 



B. Background fluxes of electron and positron 

Positrons in our galaxy are not only produced by S particle annihilations but also by the scattering of cosmic-ray 
protons with the interstellar medium |48j . The positrons produced from the later sources thus act as background for 
the positrons produced from DM annihilations. The background positron fraction can be defined as 



bkg 
sec, e+ 



bkg 



$ bkg 



bkg 



bkg 



(28) 



" prim, e~ 

where the primary and secondary electron fluxes, as well as the secondary positron flux, can be parameterized as 

O.lGe- 1 - 1 



bkg 

prim, e - 
bkg 



^GeV-'cm- W" 1 , 

15 



l + lle°- 9 + 3.2e 215 
0.70e - 7 



GeV" i cm" 2 s" 1 sr-' L , 



4> 



l + lle^ + GOO^ + SSOe 4 - 2 
4.5e - 7 



sec, c + 1 + 650 £ 2-3 + I500e 4 - 2 



GeV" 1 cm- 2 s _1 sr- 1 , 



(29) 



where e=E/(l GeV) is a dimensionless parameter. 



Results and Discussions 



The net positron flux in the galactic medium is given by 



($e+)Gal = ($ e +)bkg + * e +(^,^0)- 



(30) 



The second term in the above equation is given by Eq. ( |27[ ), which depends on various factors: b(E). \o(E : E'), 
I(Xd(E, E')), v e +, (ns)o and the injection spectrum dN e + JdE' . The energy loss (due to inverse Compton scattering 
and synchrotron radiation with Galactic magnetic field) term b(E) is determined by the photon density, and the 
strength of magnetic fields. Its value is taken to be b(E) = 10~ 16 £ 2 GeVs -1 |3S]. The number density of S DM 
in the Solar System is given by (ns)© = Pq/Ms, where pq ps 0.3GeV/cm 3 . In the energy range we are interested 
in, the value of v e + is taken approximately to be c, the velocity of light. The values of diffusion length Xd(E,E') 
and the corresponding halo function I(Xd(E, E')) are based on astrophysical assumptions [451 146|. By considering 
different heights of the Galactic plane and different DM halo profiles the results may vary slightly. In the following, 
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the Galactic plane height is taken to be less than 4 kpc. which is referred to as the "MED" model |4"51, and we 
have used the Navarro-Frenk- White (NFW) DM halo profile (SU] 



r 



P(r)=PcA-) I R I (31) 



1 



to determine the halo function I(Xd(E, E')), where r s ps 20kpc and r s» 8.5kpc. 

We use the program DARKSUSY [51] to compute electron and positron fluxes from S annihilations S — > 
^2^2j^2 • We then determine, for 100 GeV < Mg < 1 TeV, what is the maximum annihilation cross 

section allowed for the fluxes not to exceed PAMELA observations. In this range, we found the approximate empirical 
upper bound: 

/ if \ 1-95 

<S e ) < 1800 x • (32) 



The constraint Eq. (32 1 only tells us that bigger boost factor are excluded. If one wants to fully account for the 
anomalous cosmic ray fluxes through DM annihilations, the boost factor gets lower-bounded, again in the 100 GeV < 
Ms < 1 TeV range: 

" m s ^ 1M 



(S e ) > 1000 x ( ) (33) 



If both Eqs. (32 1 and (33 1 are fulfilled, then cosmic ray fluxes are within 1 standard deviation of PAMELA data. 
Of course, if all PAMELA data have to be explained, A/5 > 400 GeV is required, given the annihilation channel 
S'S —> /12/12, h-2 — > - 

Actually several constraints exist on large annihilation cross sections, relying on different physics, but all sensitive 
to DM annihilation products. When compared to the fiducial value of the annihilation cross section, these constraints 
apply in turn to the Sommerfcld enhancement (S e ). At high redshift, the energy deposition of the charged leptons 
may induce perturbations of the cosmic microwave background photon spectra [52_ , reionization and heating of the 
intergalactic medium |53| . providing strong constraints. At the recombination time during which DM relative velocity 
is j3 ~ 10~ 8 , a bound on S e is inferred |52| : 

5 ^ 480 X (l^)- (34) 

Stringent constraints arise from inverse Compton gamma rays in the Galaxy. The muons produced in DM annihilations 
subsequently decay into electrons. This population of electron yields irreducible high-energy gamma rays through 
inverse Compton on the Galactic radiation field. We consider the Fermi data released in [SI], at galactic latitude 
|b| > 10° and the analysis of Ref. [55] • For a NFW profile and a final state into 4/x, the allowed boost factor is ~ 300 
for a dark matter mass of 400 GeV. Notice that Ref. |56j considers a Galactic latitude closer to the Galactic center 
and is therefore more sensitive to the DM density profile. In this case and for a NFW profile the maximum allowed 
boost factor at Ms = 400 GeV is ~ 100. The model is thus in great tension with the Pamela anomaly. However 
with an isothermal profile instead, a boost factor up to ~ 1000 is allowed, attenuating the constraints on the model 
parameter space. As for the extragalactic gamma ray constraints, discussed in Rcfs. [55, 57 , they strongly depend on 
the assumptions on the history of structure formation. It turns out that the parameter space we are considering is 
allowed for a conservative choice of the halo concentration parameter, see Fig. 4 of [55] and Fig. 6 of Ref. [57] (which 
considers a two-muon final state case). For the cases we consider, constraints from recombination and from diffuse 
gamma rays are of the same order. In Figs.(|3|-(|5|, we depicted the boost factor satisfying the constraints Eqs.( 32p4 l 
within a green (gray) band. 

In Fig. @, we show the comparison between the positron fraction obtained from S annihilations with the positron 
fraction observed by PAMELA, AMS and HEAT, for the two typical points defined in Table|T]at M s = 400 GeV. The 
quoted values of the couplings are inferred from Fig. ^ , and provide a relic density in the WMAP7 range as well 
as saturate current direct detection bound. The first set of parameters is in a good agreement with the anomalous 
positron fraction observed by PAMELA. However, the corresponding boost factor (S e ) ~ 200 is at the border line of 
gamma ray constraints coming from FermiLAT [55j-|57 and reionization [53J. The second set cannot entirely account 
for the excess measured by PAMELA since the largest boost factor allowed by direct detection is about 8. We conclude 
that in order to fully explain the positron fraction together with an observable cross section on nucleon, bigger values 
of it x fst/>, and conversely lower values of the mixing angle are mandatory. 
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TABLE I: Typical points used for the analysis of the PAMELA positron fraction. The DM mass is fixed at 400 GeV. 
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FIG. 7: Electron plus positron flux from S*S — > 2(i + /i at Ms = 400 GeV. The red solid line (black dashed line) stands for 
the illustrative point 1 (2) quoted in Table [I] 



Concerning the electron plus positron flux observed by FermiLAT, Ms ~ 1 TeV is necessary if the flux stems only 
from DM annihilations. However, for such high masses, the required boost factor is (S e ) > 10 3 . Such a high value is 
in great tension with reionization constraints [55]. In Fig. ([7]), we compare FermiLAT, HESS, PPB-BETS and ATIC 
data on electron plus positron fluxes with the predictions of this model for the representative points mentioned in 
Table|l] for which M s = 400 GeV. The first point, depicted by a red solid line, yields (S e ) = 195 and can fit PAMELA 
data as well as the low energy flux observed by FermiLAT. Of course, with Ms = 400 GeV, not all the energy range 
can be explained. The second point, depicted by a dashed black line and for which (S e ) — 8, cannot account neither 
for the PAMELA results nor for the FermiLAT ones. In both cases, FermiLAT electron plus positron flux can only 
explained by adding an extra source of astrophysical origin. 

VI. CONCLUSIONS 

In this paper we studied a hidden Abelian extension of the standard model. The DM is a complex scalar, singlet 
under the SM gauge group but charged under the hidden sector. We also introduce a light scalar $, whose non-zero 
vev u breaks U(1)h to a Z 2 symmetry under which all fields but S are even. As a result S is a stable DM candidate, 
with mass ranging from the GeV to the TeV scale. The three fields H, $ and S couple together via three Higgs portal 
couplings. The physical scalars h and <f> mix together with the mixing angle 9 m oc /h4>- 

The relic density of S mainly results from the annihilation channel S — > h 2 h 2 through fs<j>- AH three portal 
couplings enter in the spin-independent DM-nucleon cross section. While the usual Higgs-mediated channel depends 
on Ihs- the mixing between h and 4> provides an additional channel, mediated by the light scalar /12, which is 
oc u x fs^ x 8 m . Given the mass scales we consider, the main contribution to the direct detection signal actually 
comes from this mixing term. For the parameter space we scanned over, solutions are found saturating or exceeding 
current experimental bounds from XENON10 and CDMS-II, or are in the reach of sensitivity of XENON100. The 
model also provides indirect signatures of S through cosmic ray flux measurements. The main annihilation channel 

S — > h 2 h2, followed by fast h 2 — > decays, ends up in high-energy electron and positron fluxes. The suppression 

of S annihilations into standard Higgs compared to our dominant channel, together with the light mass of h 2 , entail 
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that no antiproton flux is expected at a significant level. A Sommerfeld enhancement of the current DM annihilation 
cross section occurs through the light hi exchange, which explains cosmic ray excess observations. This enhancement 
depends on u x fg^. 

Interestingly, in this model, direct and indirect DM searches constrain same part of the parameter space. More 
precisely, in order to fully account for the anomalous positron fraction observed by PAMELA, large values of u x fg^ 
are required. Such large values give rise to large direct detection signals, saturating the current experimental exclusion 
limits. From this, an upper bound on the mixing angle is inferred. As an example, for Mg — 400 GeV and u x fg^ ~ 
410 GeV, 9 m should be less than 1.9 x 10~ 6 in order to satisfy simultaneously the current direct and indirect detection 
limits. We stress that in all models where the Sommerfeld enhancement occurs thanks to a light scalar that mixes 
with the Higgs particle, direct and indirect detection of dark matter are tightly connected. 
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